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A long-lived qubit is usually well-isolated from all other systems and the environments, and so 
is not easy to couple with measurement apparatus. It is sometimes difficult to implement reliable 
projective measurements on such a qubit. One potential solution is spin amplification with many 
ancillary qubits. Here, we propose a spin amplihcation technique, where the ancillary qubits state 
change depending on the state of the target qubit. The technique works even in the presence of 
inhomogeneous broadening. We show that fast and accurate amplification is possible even if the 
coupling and frequency of the ancillary qubits is inhomogeneous. Since our scheme is robust against 
realistic imperfections, this could provide a new mechanism for reading out a single spin that could 
not have been measured using the previous approaches. 

PACS numbers: 


Accurate quantum measurements and long coherence 
times are typical requirements for quantum information 
processing. It is essential to perform quantum measure¬ 
ments when we require quantum error correction[lj or 
need quantum feedback technology 0. 

Moreover, long coherence time is important in terms 
of retaining the quantum information during the waiting 
time for the next operations. A lot of progress has been 
made on the technique for fabricating a long-lived qubit 
with many physical systems, and the coherence time has 
been significantly improved. For example, an electron 
spin has a coherence time longer than a second Q, while 
the coherence time of a nuclear spin is longer than a 
minute [^. Even for an artificial two-level system such 
as a superconducting qubit, the coherence time is of the 
order of tens of microseconds 

It is sometimes difficult to read out the signal from a 
long-lived qubit, because long-lived qubits do not gen¬ 
erally couple well with the measurement apparatus. To 
increase the coherence time, the qubit should be isolated 
so that it has a tiny coupling strength with other exter¬ 
nal systems. On the other hand, to read out a qubit with 
high fidelity requires a strong coupling with the measure¬ 
ment apparatus, which may be contradicted by the fact 
that the qubit is isolated from other systems. 

One solution for this problem is spin amplification. If 
we couple the target qubit to read out with many ancil¬ 
lary qubits, we can use the ancillary qubits to increase 
the effective signals, even if there is little individual cou¬ 
pling between the target qubit and ancillary qubits. If the 
target qubit is in the ground state, the ancillary qubits 
remain in the ground state, while the ancillary qubits will 
be excited if the target qubit is in the excited state. There 
has been a lot of researches in this direction. A simple 
quantum circuit, which consists of controlled-NOT gates, 
is a well-known example • However, in this scheme 
each spin must be independently accessible and inter¬ 
act with the target qubit and it is difficult to increase 
the scale of such a scheme. Although there have been 


some proposals as regards a scalable spin amplification, 
they require a special co nfig uration for ancillary qubits 
with tailored couplings [uMl 2] and/or the protocol is frag¬ 
ile with respect to certain some imperfections such as the 
inhomogeneous broadening or thermal excitation of the 
ancillary qubits [13, Hi- 

In this paper, we suggest another way to amplify the 
signal of the target qubit where many ancillary qubits 
are collectively coupled with the target qubit. Since this 
is a simple experimental setup, there are many physical 
systems that have this configuration. A superconducting 
qubit coupled with a spin ensemble is one candidate for 
realizing our scheme |l5l4l8l | . A nuclear (or electron) spin 
coupled with many satellite nuclear (or electron) spins 
with different frequencies is another candidate [19j- 

In this paper, we address the problem of the inhomo¬ 
geneous broadening in particular [i^, [Hj. In reality, it is 
very difficult to fabricate identical qubits so that every 
ancillary qubit can have a different frequency and also 
have a different coupling strength with the target qubit. 
Interestingly, in our scheme, we can realize rapid and 
reliable spin amplification even under the influence of in¬ 
homogeneous broadening. The target qubit is coupled 
with the collective mode of the ancillary qubits. On the 
other hand, the collective mode of the ancillary qubits is 
coupled with the subradiant modes of the qubits, which 
induces dissipation. If we try to realize coherence trans¬ 
fer between the target spin and the collective mode, the 
existence of the subradiant modes induces decoherence, 
which means the fidelity of the coherent transfer de¬ 
creases. However, with spin amplification, coherence is 
not required because the purpose is to read out the target 
qubit, which eliminates the diagonal term of the density 
matrix of the target spin if the quantum measurement 
is successfully implemented. Moreover, we have shown 
that coupling between the collective mode and subradi¬ 
ant mode only induces an energy transfer from the for¬ 
mer to the latter that can be also measured as well as 
the former. Our technique can provide a reliable readout 
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of a qubit coupled with many inhomogeneous ancillary 
qubits. 

The Hamiltonian of our hybrid system composed of the 
target qubit and ancillary qubits can be described by 



FIG. 1: Schematic of our hybrid system. The central target 
qubit is collectively coupled with the ensemble of ancillary 
qubits. The target qubit can be readout via the ancillary 
qubits after the amprification technique 
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where we set h = 1. Here, az denotes the Pauli z op¬ 
erator for the target qubit and d^^+Menotes the ladder 
operator, wt and ujj denote the resonant frequencies of 
the target qubit and the j th spin, respectively, tf^ 
and dzj denote the ladder operator and Pauli z operator 
of the j th ancillary qubit, respectively. It is worth men¬ 
tioning that, as long as most of the ancillary qubits are in 
the ground state and only a small portion of the ancillary 
qubits are excited for the amplification, we can consider 
the ancillary qubit to be a harmonic oscillator. So we 
can replace a - with the annihilation(creation) oper¬ 
ator ■ Sj is a real number and denotes the coupling 
strength between the target qubit and the j th ancillary 
qubit. Xd denotes the Rabi frequency of the microwave 
applied to the target qubit. We assume that the cou¬ 
pling strength gj is much smaller than the inhomogeneous 
broadening width of this system. By eniployin^a rotat¬ 
ing frame defined by U{t) = exp[—iujdt{^f- p^j)] 

and performing a rotating wave approximation, we have 

N 
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We will now model the effect of inhomogeneous broad¬ 
ening with a Lorentzian distribution of the resonant fre¬ 


quencies of the spin ensemble as 


/(w) 


1 7/2 
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where, 7 denotes the width of the Lorentzian, and w de¬ 
notes the average spin ensemble resonant frequency. We 
define the collective mode as follows 



(4) 


where G = pYhPi 9^- Since we are mainly interested 
in the target qubit and collective mode, we can trace 
out the sub-radiant modes. In this case, the collective 
Hamiltonian can be written as 221 


H = He + Hd ( 5 ) 

He = 2 4" + ^4(7“) -I- (w — uJd)A'^A 

( 6 ) 

= Y(d+-bo--) (7) 


The width of the Lorentzian 7 can be incorporated in 
the dynamics as the decay term of the master equation, 
which we will discuss later. 

In our scheme, we will drive the target qubit with a mi¬ 
crowave pulse for the spin amplification. The collective 
mode coupled with the target qubit may differ from the 
collective mode coupled with the microwave line, which 
is called a mode mismatch[^. This problem means that 
driving the collective mode by the coupling between the 
microwave line and ancillary qubits is not straightfor¬ 
ward. On the other hand, as we explain later, it is possi¬ 
ble to drive the collective mode of the ancillary qubits by 
driving the target qubit. We have developed an efficient 
technique for calculating the number of excitations in the 
spin ensemble after driving. If we include the all freedom 
of N ancillary qubits, we need to deal with 2^+^ sized 
matrix, which is intractable. On the other hand, with 
the technique that we propose here, we need to solve a 
Hamiltonian that includes a single qubit and a harmonic 
oscillator, which is much easier to simulate. Using Eql2] 
as a basis, we can write the Heisenberg equation, 



-igjait) - i{ujj -ujd- 


( 8 ) 


where we add the energy decay term T of the ancillary 
qubit. However, this decay rate is usually negligible com¬ 
pared with that of other noise. For example, the energy 
relaxation time of an electron such as the NV” center 
exceeds 200 seconds [2J . So we adopt a limit of L —>• 0, 
later. Now 


dt 
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Now, taking the Laplace transforms of both sides of EqlHl 
under the initial conditions aj ( 0 ) = 0 , cr~ ( 0 ) = 0 , we have 




igjP{s) 


+ i{ujj -LOd- i^) 


( 10 ) 


Here, aj{s) = (t)],/3(s) = C[cj (t)]. Using the 

Lorentzian distribution m and adopting residual inte¬ 
gration, we obtain 


N 


G 


= - 


iCPis) 


1=1 


s + i{uj - uJd - 


( 11 ) 


Taking the inverse Laplace transformation of both sides, 
we obtain 

iGa~{t) = - ii^-uJd - (12) 

By using EqU^] and EqlHl in the limit of 7 >> L we 
obtain 

, IV 

— (-il(t)i(t)+ ^a](t)a(t)j) =jA^t)A{t) (13) 
1=1 

This means that the rate of increase of the subradi¬ 
ant modes (—Hl(t)H(l) -I- ®](0®l(^)) depends on 

the collective mode (Hl(t)yl(t)) and the inhomogeneous 
width ( 7 ). By integrating EqlTSl we obtain 


V a(t)]a(t)j = i(t)li(t) + 7 / A{t'AA{t’)dt' (14) 

,=i -lo 


By calculating he dynamics of A(t), we can estimate the 
excitation number X)j=i ®l( 0 ^®l (0 

We intuitively explain how our amplihcation works. 
If the detuning between the energy of the target qubit 
and the frequency of the collective mode is much larger 
than the coupling strength, we can derive the dispersive 
Hamiltonian as follows 12^ . 


the driving caused by the detuning. Throughout this 
paper, we fix the driving frequency as a; -I- 

Although we later solve the master equation numeri¬ 
cally, it is also possible to describe the dynamics of the 
collective mode in an analytical way by using an approx¬ 
imation. By diagonalizing the He in Eql 6 l we obtain 


He = Y.Hn 

n 

Hn = f 2 +(n) |-|-,n) (+,n\ + f 2 _(n) |-,n) (-,n| 


r> ^ \ t i ^ _i_ \/A^ -I- 4G2(n -I- 1) 

n±{n) = {n+ -){uj - UJd) ±- - - 


(17) 

(18) 

(19) 

( 20 ) 


The eigenstate of the Hamiltonian is 


|-|-, n) = cos((/>„/ 2 ) |e, n) + sin(^„/ 2 ) \g, n -I- 1 ) 
|-, n) = -sm{(j)n/2) |e, n) + cos{4>n/2) \g, n+\) 


(pn = tan ^ 


2GVnH^\ 

A ) 


( 21 ) 

( 22 ) 

(23) 


where |e) {\g)) denotes the excited (ground) state of the 
target qubit and \n) denotes the number state of the col¬ 
lective mode. Since we can diagonalize the Hamiltonian 
He, we can efficiently use an interaction picture where 
we define (j = and Hi{t) = HdU■ By using the 

condtion A ^ C?, we obtain 


Hi{t) Xeff(\e){e\{A + it) - |g)(g|(ie*^‘ + 

(24) 


where Ae// = denotes the effective Rabi frequency 
of the collective mode. If the initial state of the target 
qubit is |e) and it remains unchanged during the time 
evolution, we can derive the Hamiltonian for the collec¬ 
tive mode Hj'^’^^\t) ~ Ae//(i-l-it). On the other hand, 
if the target qubit is kept in the ground state, we obtain 
Hi{t) ~ —Ae//(ie®^‘ -I-ite“®^‘). By returning to 
the Schrodinger picture, we obtain 


UJT-UJd. G2 ^ - 

Hdisp = -^-cr^-I-— cr^cr +{uJ-ujd + -^(JAA'A 

(15) 


where Ad = 0 and A = [ujt — ujd) — {uj — uJd)- This means 
that, depending on the state of the target qubit, there is 
a dispersive energy shift as follows. 

G^ 

w => w ± (16) 

So, if the driving frequency is w-l-^ and the target qubit 
is in an excited state, the ensemble of the ancillary qubits 
is efficiently excited. On the other hand, if the driving 
frequency is w -I- ^ and the target qubit is in the ground 
state, the spin ensemble is not significantly affected by 


T^arac = (w — Wd ± -^)A’f A ± Ae//(A-I-A’^) (25) 

where the sign ± depends on the state of the target qubit. 
With this Hamiltonian, we obtain the Heisenberg equa¬ 
tion of A as follows, 

d '' 'y ^ 

—A = -i{u!-uJd±-^-ipA-iXeff (26) 

where we added the energy relaxation term —il^A. When 
the target qubit is in an excited state, we solve this equa- 
tion analytically with uJd = uj and obtain 

4A2 / 

(At(t)A(t)) =-e-^‘)2 (27) 











On the other hand, when the target spin is in the ground 
state, we obtain 


{A\t)A{t)) = 


^eff 




1 —2cOs( 2 *_|_g T'* 

(28) 


It should be noted that, if the effective Rabi frequency 
is far smaller than the dispersive energy shift such as 
Xeff the excitation of the collective mode is sig¬ 

nificantly suppressed when the target qubit is prepared in 
the ground state. So we retain this condition throughout 
this paper. 
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FIG. 2: Collective mode excitation number against time. 
Solid lines show analytical solutions and dashed lines show 
numerical calculations where we set 7/27r = 12.5 MHz, 
G/27r=75 MHz, A/27r = 412.5 MHz, and Ad/27r = 40 MHz. 
These are typical parameters of the hybrid system composed 
of superconducting flux qubit and NV- centers[T^. When ob¬ 
taining analytical solutions, we use a dispersive Hamiltonian 
where the state of the target qubit is not changed, which is 
a valid approximation for a short time region (dispersive re¬ 
gion). However, for a numerical simulation, the excited state 
of the target qubit decays into the ground state due to the flip 
flop term of the JC Hamiltonian and subradient mode, which 
leads a difference from the analytical approach for a longer 
time scale. 

We calculate the time evolution of the collective mode 
by solving the following master equation. 

^ = -^[H,p] + ^[2A0 - A^Ap - 0A] (29) 

It is worth mentioning that the inhomogeneous broad¬ 
ening of the frequency can be considered an energy loss 
of the collective mode [1^ . In Fig [21 we plot the time 
evolution of the collective mode by using the analytical 
approach described above and a numerical approach to 
solve the master equation of Eq. [29l We assume that the 
initial state is |e,0) or |g, 0). In the analytical approach, 
the state of the target spin does not change during the 
time evolution, because we use a dispersive Hamiltonian. 


FIG. 3: Spin amplification process in our scheme. Depending 
on the initial state of the target qubit (|e) or Ig)), the ancillary 
qnbit excitations numbers differ, and we plot this difference 
against time. We set the parameters at G/27r = 75 MHz, 
Xd/2n = 40 MHz, A/2tv = 412.5 MHz 


When the initial state of the target qubit is an excited 
state, the collective mode of the spin ensemble is excited 
efficiently in both the numerical simulation and the ana¬ 
lytical calculation for a short time scale. However, since 
the energy of the target qubit is gradually transfered to 
the subradiant modes of the ancillary qubits, the target 
qubit decays into the ground state and so the collective 
mode will not be efficiently excited over a long time pe¬ 
riod in the numerical simulation . This induces a 
disagreement between the analytical solution and numer¬ 
ical simulation in this regime. On the other hand, when 
the initial state of the target qubit is the ground state, 
the state of the target will not be significantly changed, 
and so there is an excellent agreement between all t val¬ 
ues of the analytical solution and numerical simulation 
where the excitation of the collective modes is efficiently 
suppressed. 

We plot the difference in the excitation population of 
the ensemble when the target qubit is in its excited and 
ground states in Fig [3] Although the inhomogeneous 
broadening reduces the efficiency of the amplification, we 
can still obtain a difference of more than 10 excitations 
in the ancillary qubits for 7/27r = 10 MHz. This means 
that we will have a 10 times better signal with our ampli¬ 
fication technique than without any amplification, which 
would provide us with a practical tool with which to read 
out a single qubit such as an electron spin or a nuclear 
spin. 

In conclusion, we have proposed a spin amplification 
technique that can be used in the presence of the inhomo¬ 
geneous broadening. We considered a simple experimen¬ 
tal system where a target qubit is collectively coupled 
with many ancillary qubits. We showed that, even if the 
ancillary qubits are inhomogeneous, we can amplify the 
signal from the target spin, and obtain a gain of 10 in 
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such a system by applying a continuous microwave. Our 
technique opens the way to the fast and robust readout 
of a single spin. 
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